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Dense quantum plasmas are ubiquitous in planetary interiors and in compact astrophysical objects 
(e.g. the interior of white dwarf stars, in magnetars etc.), in semiconductors and micro-mechanical 
systems, as well as in the next generation intense laser-solid density plasma interaction experiments 
and in quantum x-ray free-electron lasers. In contrast to classical plasmas, one encounters extremely 
high plasma number density and low temperature in quantum plasmas. The latter are composed of 
electrons, positrons and holes, which are degenerate. Positrons (holes) have the same (slightly differ- 
ent) mass as electrons, but opposite charge. The degenerate charged particles (electrons, positrons, 
holes) follow the Fermi-Dirac statistics. In quantum plasmas, there are new forces associated with 
i) quantum statistical electron and positron pressures, ii) electron and positron tunneling through 
the Bohm potential, and iii) electron and positron angular momentum spin. Inclusion of these 
quantum forces provides possibility of very high-frequency dispersive electrostatic and electromag- 
netic waves (e.g. in the hard x-ray and gamma rays regimes) having extremely short wavelengths. 
In this review paper, we present theoretical backgrounds for some important nonlinear aspects of 
wave-wave and wave-electron interactions in dense quantum plasmas. Specifically, we shall focus on 
nonlinear electrostatic electron and ion plasma waves, novel aspects of 3D quantum electron fluid 
turbulence, as well as nonlinear ly coupled intense electromagnetic waves and localized plasma wave 
structures. Also discussed are the phase space kinetic structures and mechanisms that can generate 
quasi-stationary magnetic fields in dense quantum plasmas. The influence of the external magnetic 
field and the electron angular momentum spin on the electromagnetic wave dynamics is discussed. 
Finally, future perspectives of the nonlinear quantum plasma physics are highlighted. 

PACS numbers: 05.30.Fk,52.35.Mw,52.35.Ra,52.35Sb 
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I. INTRODUCTION 



The field of quantum plasma physics has a long and diverse tradition [H [21 |3j HI |5] , and is becoming of increasing 
current interest [H |7], motivated by its potential applications in modern technology (e.g. metallic and semiconductor 
nanostructures-such as metallic nanoparticles, metal clusters, thin metal films, spintronics, nanotubes, quantum well 
and quantum dots, nano-plasmonic devices, quantum x-ray free-electron lasers, etc.). Due to the recent development 
of ultrafast spectroscopy techniques, it is now possible to monitor the femtosecond dynamics of an electron gas confined 
in metallic plasmas. In dense quantum plasmas, the number densities of degenerate electrons and/or positrons are 
extremely high, and the plasma particles (mainly electrons and positrons) obey Fermi-Dirac statistics. 

The quantum degeneracy effects start playing a significant role when the de Broglie thermal wavelength Xb = 
h/{2iTmf,kBTy^'^ for electrons is similar to or larger than the average inter-electron distance Ue i.e. when [6l|7] 



n, 



■M (1) 



or, equivalently, the temperature T is comparable or lower than the electron Fermi temperature Tpe = Ep/kB, where 
the electron Fermi energy is 

Ep = ^i3n')V^nl/^ (2) 
Znip 



so that 



= i(37r2)V3(n,A|)V3 > 1. (3) 



T 2 

Here h is the Planck constant divided by 27r, rie is the electron number density, me is the rest electron mass, and fc^ 
is the Boltzmann constant. 

When the temperature approaches the electron Fermi temperature Tpe, one can show, by using the density 
matrix formalism |H], that the equilibrium electron distribution function changes from the Maxwell-Boltzmann 
oc exp(— i^/fc^T) to the Fermi-Dirac distribution cx {2/h^)[exp{{E + ii)/kBTFe) + ^~ , where E is the electron 
kinetic energy and is the chemical potential. In a dense Fermi plasma, the Thomas-Fermi screening radius reads 

A. = ^ (4) 

which is the quantum analogue of the Debyc-Hiickel radius. Here the electron Fermi speed 

Vp, = {2EF/m,y^' ^ —{STr'n.y/^ (5) 
me 

is the speed of an electron at the Fermi surface. 

A measure of the importance of collisions in a dense plasma is the quantum coupling parameter, which is the 
ratio between the interaction energy Eint = e^rifJ^ and the average kinetic energy E'^m of electrons, where e is the 
magnitude of the electron charge. For a classical plasma, the kinetic energy is fcsT and hence we have Fc = Eint/kBT 
in the classical case. In a quantum plasma, we have instead Ekin = Ep, which gives the quantum coupling parameter 

[siaii] 

Ep (3^2)2/3^2^1/3 ^ |^„^A3,; ^y2kBTpJ-''' 

(where we have left out proportionality constants for the sake of clarity) is analogous to the classical one when 
Xp ^ Xp)- The various plasma regimes are illustrated in Fig. n] where the straight lines correspond to i) the 
distinction between the classical and quantum plasmas, x = 1; ii) the limit between coUisionless and coUisional 
classical plasmas, Tc = 1, and iii) the limit between coUisionless and coUisional quantum plasmas, Fg = 1. In Fig. 
[l] various experimental and naturally occurring plasmas are exemplified. Also indicated in Fig. [l] are the kinetic 
equations used to model the plasma in each regime. The Vlasov and Wigncr equations are thus used to model 
coUisionless plasma in the classical and quantum limits, respectively, while "Boltzmann" indicates coUisional kinetic 
models in a classical plasma. The kinetic models for a coUisional quantum plasma is labeled "Wigner (-fcoU)" in 
Fig. [l] We note that a quantum plasma becomes coUisionless Fq < 1 when the mean distance between electrons is 
of the order of the Bohr radius oq, i.e. d = < [{Stt'^)'^^^ /8'!r]aQ w 0.38ao, where oq = /rriee^ « 0.53 A. This 
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corresponds to a number density of the order ^ 1.22 x 10^^ cm""^, which is three orders of magnitude larger than 
the electron density in a typical metal, for example gold has Ue = 5.9 x 10^^ cm~'^ at room temperature. 

In Fermi-degenerate matters, there is an effect called the Pauli-blocking, which strongly reduces the electron-electron 
and electron-ion coUisional rates. Namely, at moderate temperatures, only electrons within an energy shell of thickness 
ksT about the Fermi surface (where the electron energy equals Ep) can undergo collisions. For these electrons, the 
electron-electron collision rate is of the order ksT /Ti, and the average collision rate is obtained by multiplying this 
expression by T/Tp. The resulting collision frequency v^,^ divided by the electron plasma frequency is |7] 



Ep 

huJri 



T y 



TpJ Y^/"^ \TfJ 



T y 



(7) 



Hence, v^.^ <C ojpe when T <^ Tp and Tq > 1, which is relevant for metallic electrons. For example, at room 



temperature we have Ve 
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interactions, 



10 



16 < 



which is much smaller than the typical coUisionless frequency of collective 



Also, the typical collision rate for electron-lattice (ion) collisions. 
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smaller than iVpg by one order of magnitude. Therefore, a coUisionless regime seems to be relevant for free electrons 
in a metal on a time-scale of the order of a femtosecond. In denser plasmas, such as in stellar interiors [101 EJ [12] 
and in incrtial fusion schemes [T31 HH [15] , the relative effects of collisions decrease even further and leads to increased 
electron transport and heat conductivity. 
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FIG. 1: Schematic plasma diagram in the logjQ(r) - log(n) plane: lONO: ionospheric plasma, SPACE: space plasma, CORONA: 
solar corona, DISCHA: typical electric discharge, TOK: tokamak/magnetic fusion experiments, ICF: inertial confinement fusion, 
MET: metals and metal clusters, JUP: Jupiter's core, DWARF: white dwarf star. After Refs. [GUT]. 

More than sixty five years ago, Wigner |T6l introduced a phase-space formalism to treat a quantum state of charged 
particles in a coUisionless quantum system. He introduced a quantum distribution function of the phase-space variables 
f{x,p.t). The Wigner distribution function is defined as 



f{x,p,t) = 



(2^) 



3Af 



T,x + — r, t I exp{—ipT)dT, 



(8) 



where N is the number of particles in the system, x{x;1,X2---,xn} said p{pi,p2, Pn} are the sets of coordinates 
and momenta for the particles, t is the time, and p{x,x',t) is the density matrix. The Wigner function is not 
a probability density in phase space x,p, because it can take negative values. The Wigner distribution function 
Wi{x,p.t), corresponding to the wave function ipi(x,t), can be expressed as 



W^{x,p,t) 



dy (V'j(a; - y/2, t)ip*{x + y/2, t)) exp{-py/h,), 



(9) 
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which has the property 



dpW^{x,p.t) =< \iJ^{x,t)\'^ >, 



(10) 



where the asterisk denotes complex conjugate. The quantum kinetic equation based on the Wigner distribution was 
developed by Moyal [T7], and is now referred to as the Wigner-Moyal description for treating statistical effects on 
electron plasma waves in a quantum plasma [18, . 

Analytical investigations of collective interactions between an ensemble of degenerate electrons in a dense quantum 
plasma dates back to early fifties. Specifically, KHmontovich and Silin [T] and Bohm and Pines [2J [31 IH IS] presented 
properties of linear electron plasma oscillations (EPOs) in a dense quantum plasma. In the latter, electrons, positrons, 
and holes are degenerate, while ions are cold and classical (typically, in a dense quantum plasma the ion Fermi speed 
is much smaller than the electron Fermi speed). Accordingly, electrons, positrons and holes have a Fermi-Dirac 
distribution function [19' , contrary to the Boltzmann-Maxwell distribution function for charged particles in a classical 
plasma. 

The dispersion relation for high-frequency electron plasma oscillations in a dense quantum plasma with fixed ion 
background reads (see Appendix B) 



1 



47re^ 



/o(u) 



Ami 



d^u = 0, 



(11) 



which was also obtained by Bohm and Pines [3] by performing a series of canonical transformations on the Hamiltonian 
of the system of individual electrons, interacting via the electrostatic force. In the zero temperature limit, we have 
(see Appendix B) 



pe 



fikV, 



Fe 



V, 



Fe 



flkVpe 



fik 
2™. 



log 



hk 



log 



hk 
2m ^ 

hk 
2m, 



VFe 



hk 
2m, 



hk 
2m^ 



(12) 



= 1 + Xe = 0, 



which, in the limit Tik/nie — > yields 



2kVF 



■log 



kVpe 



kVp 



= 0, 



(13) 



where we have assumed that uj is real and uj/k > Vpe- Here u) is the wave frequency, k is the wave vector, and 
(47rnee^/me)^/^ is the elec tron plasma frequency. On the other hand, for small wavenumbers up to terms 



containing fc , we obtain from ( 12 1 



:k^V, 



(1 



(14) 



4m^ ' 

where a = (48/175)m^ V^^/Ti^w^^ w 2.000(aj]no)^^^, where gq = Ti^ /niee^ w 53 x 10"^° cm is the Bohr radius. 
Equation (14) shows that the wave dispersion arises due to the finite width of the electron wave function in a dense 
Fermi plasma [H |20l ETJ |22l [23l [24] . 

Furthermore, in the low phase speed limit, viz. lo <^ kVpe, the dielectric constant for ion oscillations reads 



e(^,k)«l 



3r A;4/4to2 



(15) 



which yields, with e(a;, k) — 0, the ion oscillation frequency 



UJ 



(l + Q)i/2' 



(16a) 



■^pi 



{me/miY/'^LOpe is the ion plasma frequency, and rrii is the ion mass, Q = 3ujp^/{k^Vp^ + 3h^ k'^ / Ami) 



where 

and a ^ 1. For Q ^ 1, we have from (16b) 



UJ : 



kCps 1 



2f,4 



n^k 

Arriiu 



pe 



1/2 



(16b) 
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where Cps = (TVe/Smi)-'^/^ is the sound speed. 

Dispersion properties of electrostatic waves in an unmagnetized dense quantum plasmas with arbitrary electron 
degeneracy have been presented by Maafa [53] and Melrose [25]. By using the random phase approximation, the 
permittivity of a degenerate coUisionless plasma is given in textbooks [571 HH]- Furthermore, since the pure elec- 
tromagnetic wave in a non-streaming unmagnetized dense plasma do not accompany density fluctuations, the wave 
frequency is a; = {ijJ^^ + fc^c^)^/^, where c is the speed of light in vacuum. Theoretical studies of quantum statistical 
properties of dense plasmas in the presence of electromagnetic waves appear in Kremp et al. [29| and in textbooks 
[301 131j , while quantum parameter regimes are discussed by Bonitz [6 . In a magnetized dense quantum plasma, one 
finds that the external magnetic field significantly affects the dynamics of degenerate electrons and positrons, and 
subsequently there appear new collective phenomena associated with the electron angular momentum spin [32j [33] , 
the electron spin magnetic moment |34[ , and quantized Landau energy levels |35j of the Fermions in a strong magnetic 
field. It turns out that the thermodynamics and kinetics [33], as well as the dispersion properties of both electrostatic 
and electromagnetic waves [37] |3Hl [Ml HOI HH HH IH] in a quantum magnetoplasma are significantly different from 
those in an unmagnetized quantum plasma. 

It was early recognized that the underlying physics of nonlinear quantum-like equations can be better understood 
by casting those equations in the form of hydrodynamical (or the Euler) equations, which essentially represent the 
evolution of quantum particle densities and momenta. This was elegantly done by Madelung |44j and Bohm [2, by 
introducing an eikonal representation for the wave function evolution in the non-stationary Schrodinger equation. 
The derivation of the Madelung quantum fluid equations for the Pauli equation with the quantum particle angular 
momentum spin was presented by Takabayashi [iSKinillT], Bohm et al. [151112] ? Janossy and Ziegler-Naray [SD], and 
others [51|. To incorporate relativistic effects into the quantum fluid formalism, Takabayashi derived the quantum 
electron fluid equations for the Klein-Gordon equation [55] and for the Dirac equation 53, 54, 55 . Extensions have also 
been done to fluid descriptions of the Weyl equation for massless spin- 1/2 particles (neutrinos) by Bialynicki-Birula 

Recently, there has been growing and vibrant interests in investigating new aspects of quantum plasma physics 
by developing non- relativistic quantum hydrodynamical (QHD) equations [7) [SD] [^T] [S7]. The latter include the 
quantum statistical electron pressure and the quantum force involving tunneling of degenerate electrons through the 
Bohm potential [2U]. The Wigner-Poisson (WP) model has also been used to derive a set of non-relativistic quantum 
hydrodynamical (QHD) equations [7] [21] for a dense electron plasma, assuming immobile ions. The QHD equations 
are composed of the electron continuity, electron momentum and Poisson equations. The quantum force [7, 20, 21] 
appears in the non-relativistic electron momentum equation through the pressure term, which requires knowledge 
of the Wigner distribution for a quantum mixture of electron wave functions, each characterized by an occupation 
probability. Quantum transport models similar to the QHD plasma model have also been used in superfluidity |58j 
and superconductivity |59j . as well as in the study of metal clusters and nanoparticles, where they are referred to as 
nonstationary Thomas- Fermi models j60j . 

The electrostatic QHD equations are useful for studying collective interactions (e.g. different types of waves, 
instabilities, quantum fluid turbulence and nonlinear structures [Tj [21] [22l [BT] [62l [63] [64] [65] [Ml |67l (BS", ITO in dense 
quantum plasmas. The quantum kinetic and QHD equations have also been generalized to include the electromagnetic, 
ambient magnetic field and electron angular momentum spin effects [ID] [70] [7TJ [75] [73] [71] . The latter give rise to 
high-frequency spin waves, which can be excited by neutrino beams in supernovae |70|, 175] . Furthermore, studies of 
numerous collective interactions in dense plasmas are relevant in the context of i) intense laser-solid density plasma 
experiments [Tl[Tl[7SJ[77][75][7nj[Hni[HI][Hll[H31[Hl[HS], where one would be exploring new frontiers in high-energy 
density physics [SB]; ii) in the cores of giant planets and the crusts of old stars [ST] [SS] [H] ; iii) superdense astrophysical 
objects |90] [91] [92] [93] [94] [95] [96] [97] [98] (e.g. interiors of white dwarfs and magnetospheres of neutron stars and 
magnetars); iv) as well as for micro and nano-scale objects (e.g. quantum diodes [W, 1001 1101111021 [103. 104 , quantum 
dots and nanowires |105j . nano-photonics |106II107| . plasmonics [108j, ultra-small electronic devices [109. ,11011111] . and 
metallic nanostructures |112j ): v) micro-plasmas |113j . and quantum x-ray free-electron lasers |114LI115] . Furthermore, 
it should be stressed that a Fermi degenerate dense plasma may also arise when a pellet of hydrogen is compressed 
to many times the solid density in the fast ignition scenario for inertial confinement fusion |15l 11161 I117j . Since 
there is an impressive developments in the field of short pulse petawatt laser technology, it is highly likely that such 
plasma conditions can be achieved by intense laser pulse compression using powerful x-ray pulses. Here ultrafast 
x-ray Thomson scattering techniques can be used to measure the features of laser enhanced plasma lines, which will, 
in turn, give invaluable informations regarding the equation of state of shock compressed dense matters. Recently, 
spectrally resolved x-ray scattering measurements [55] [SI] have been performed in dense plasmas allowing accurate 
measurements of the electron velocity distribution function, temperature, ionization state, and of plasmons in the 
warm dense matter regime jllSj . This novel technique promises to access the degenerate, the closely coupled, and 
the ideal plasma regime, making it possible to investigate extremely dense states of matter, such as the inertial 
confinement fusion fuel during compression, reaching super-solid densities. 
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In this review article, we present the theoretical progress that has been recently made in the area of collective 
nonlinear interactions in coUisionless dense quantum plasmas. The manuscript is organized in the following fashion. 
In section 2, we shall briefly recapitulate the hydrodynamic representation of some quantum-like models that appear 
in different branches of physics. The governing equations for nonlinearly interacting electrostatic waves in an unmag- 
nctized quantum plasma are derived in section 3. Section 4 presents numerical studies of nonlinear electron plasma 
wave excitations in the form of quantized one-dimensional dark solitons and quantized two-dimensional vortices. The 
model used here is the nonlinear Schrodinger equation for the dispersive EPOs, coupled with the Poisson equation 
for the electrostatic potential. This model is also used for studying 3D quantum electron fluid turbulence in section 
|Vj where we find non-Kolmogorov-type turbulence spectra. In section VI we present recent results concerning the 



phase space (kinetic) turbulence, by using the Wigner and Vlasov models for the electron distribution function. A 
theoretical model for the generation of quasi-stationary magnetic fields in a dense quantum plasma due to the Weibel 
instability is presented in section |VII A[ The magnetization of a dense plasma in the presence of a large amplitude 



electromagnetic wave is demonstrated in subsection VII B The dynamics of electromagnetic waves in a dense magne- 



toplasma are discussed in section [VIII| Here we focus on spin waves propagating across the magnetic field direction, 
and develop nonlinear equations for low-phase speed (in comparison with the speed of light) electromagnetic waves in 
a dense quantum magnetoplasma. Finally, section [Tx] highlights our main results and describes the future prospectives 
of the quantum plasma physics research. 

II. FLUID REPRESENTATION OF QUANTUM-LIKE MODELS 

This section is included to show how different types of quantum-like models can be cast in the form of hydrodynamic 
equations. 

First, we consider the non-stationary nonlinear Schrodinger equation(NLSE) 

ot 2m 

where "ipiv^t) is the macroscopic wave function, m is the particle mass, and ?7o(|V'P) is an effective potential. The 
NLSE also arises in various physical context in the description of amplitude modulated nonlinear waves in fluids 
[1 19] 1120] ■ in transmission lines |121j . in nonlinear optics for ultra- fast communications |122[ 11231 1124] . in plasmas 
PSI \Ym Tm [T2HJ Hi inn] , and in many other areas of physics jlSJ HSl [133] . 
Introduce the Madelung transformation [ll] 

^(r,t) = ^/^exp(z^^), (18) 

where n and ipq are real, and obtain from (17) a pair of quantum hydrodynamic equations composed of the continuity 
and momentum equations, respectively, 

+ V-(nv)=0, (19) 



and 

d 



m 



dt 



vV v = -V[t/o(n) + i7s]. (20) 



dt 

Here n = n(r, t) = corresponds to the local density per unit length, and ?iV(/Jg(r, t) — rrw. The quantum potential 
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-2 



Ub = —^ (21) 

Zra \Jn 

We note that the quantum particle number density n and the quantum velocity field v can be written as, respectively, 

n(r,t) = = IV'p, (22) 

and 



h. {tlj*Vtp - V'VV'*) _ ifi ^ 
2im 2m 



In 



(23) 
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where the asterisk denotes the complex conjugate. The quantum-Hke velocity, given by (23), is a potential field, 
namely. 



V X V = (24) 

everywhere in a single-connected region. 

Let us now define the generalized vorticity on the weighted velocity field as |68j 

Vx(|V;pv) , Vl^pxv 

where the first term in the right-hand side in (25) represents the ordinary vorticity. It is well known |51l 1134] that 
in the condensate state all rotational fiow is carried by quantized vortices (the circulation of the velocity around the 
core of each such vortex is quantized). In the absence of quantized vortices, the first term is zero in view of (24). In 
such a situation, the second term in (25) determines the generalized vorticity. Various aspects of quantized vortex 
dynamics and superfiuid turbulence appear in Barenghi et al [134j. Bewley et al |135j have presented a technique for 
visualization of quantized vortices in liquid helium. 

Second, the nonlinear Schrodinger equation can be generalized by including the trapping potential 



T// \ ( 2 2, 22, 2 2\ 



(26) 



which confines identical bosons ^136j in the harmonic trap of an ultracold quantum system. Here rrib is the boson 
mass, and lJx, and lUz are the harmonic frequencies of the bosons along the x,y, and z directions, respectively, of 
a Cartesian coordinate system. The nonlinear dynamics of the Bose- Einstein condensates (BECs) |136l I137j is then 
governed by the Gross-Pitaevskii equation [138. ,139. ,140) 



ih 



d(l){T, t) 



dt 



2m 



VV(r,t) - 14V(r,i) - G|VP(r,i)V(r,i) = 0, 



(27) 



where for the BECs we have L'odV'P) = (47r?i^as/m;,)^/^|'!/;p = Gji/'P- Here mf, is the mass of the bosons and as is 
the scattering length for boson-boson collisions. The BECs are repulsive (attractive) for G > (<)0. 

Introducing the Ansatz ii){r,t) — ^ ^^(r, t) exp[</5b(r, t)] in (27), we obtain generalized quantum hydrodynamic 
equations [T4T|[T42] 



~dt 



V • (nuh) = 0, 



and 



where the particle fiux 



dt 



-V V, 



2mby/nb 



nb{r,t)ub{r,t) 



2inib 



Gub 



(7/;*V?/' - VVV'*) 



(28) 



(29) 



(30) 



with Uf, = (?i/mf,)V(^(,(r, t). Equation (29) establishes the irrotational nature of the superfiuid motion of the BECs. 
Equations (29) and (30) can be used to study the linear and nonlinear properties of BECs. 

Next, we consider the dynamics of a non-relativistic single Fermi (1/2— spin) particle (a degenerate electron) 
governed by the Pauli equation |143l I144j 



ih 



5* 



dt 2m. 



2mpC 



{A-y + v -A) 



2meC^ 



^ecr ■ B 



* = 0, 



(31) 



where \E'(r, <,cr) is the wave function of the single particle species having the spin s = l/2cr, <t is the Pauli spin 
matrices, A is the vector potential, is the scalar potential, B = V x A, and fj,B = eh/2meC is the Bohr magneton. 
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By using the Madelung representation for the complex wave function |145j 



^'(r,t,cr) = s^yne{r,t,a■)exp 



iSe{r, t, a) 



one can obtain from (31) the quantum magnetohydrodynamic equation |19l 1145] 

'npe 



due 
dt 



V • 



= 0, 



(32) 



(33) 



and 



d 



1 



dt m, 



-Pe • V p 



IdA 1 

^<^+ "IxT - X V X A 
c at c 



^(Y:^\^^,^^^^.^)^ (34) 



2m 



where s mimics the spinor through which the electron-1/2 spin properties are mediated, ne(r, i, cr) = vj/v]/* represents 
the probabihty density of finding the single electron at some point with a spin s. We have denoted the generalized 
electron momentum = V^e — ihs*\J s + (e/c)A. 

We can now express the quantum electron velocity [145] 



= ^ + i s*Vs ■ 



2top 



the spin density vector 



s — — s crs. 
2 



and the spin vector transport equation [45j 



ds e , „, 1 
sxB + 

at rUpC iTLpUp 



ds 

s X - — I Ue-— 
dxi, \ dxk 



(35) 



(36) 



(37) 



where we have used the summation convention for repeated indices, and have denoted d/dt = (d/dt) + Ve • V is the 
total derivative. The electron momentum and electromagnetic fields are coupled via the Maxwell equations. 

Ignoring the electromagnetic fields and the particle spin, one can obtain, after linearizing (33) and (34), the frequency 
of electron oscillations 



~ 2m, ' 



(38) 



where k is the wave number. 

Finally, we consider interaction of an electron with both background electrons and singly charged positive ions. 
The electron dynamics is governed by |146j 



(39) 



where 



UM')^e'l d3r'_L_(|^(r',i))|2_„^(r',i)), (40) 

is the potential responsible for the interaction of an electron with background matter which includes electrons and 
positively charged ions with the number density ni(r',t). The wave function is normalized on the number density 
of electrons, viz. ne{r,t) — \ip{r,t)\'^. We note that Eq. (39) accounts only for the Coulomb interactions between 
electrons and ions, and completely ignores the quantum statistical pressure, the self-consistent ambipolar field arising 
from the charge separation, and the electron spin- 1/2 effect. 
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Assume that tp{r,t) = ipo + ipi{r,t), where \ipo\'^ = no represent the unperturbed electron number density and the 
perturbation wave function ipi(r,t) for spherically symmetric oscillations has the form (Afe/r)sin(kr) exp(— iwt), with 
Ale being the normalization constant. Thus, the dispersion relation deduced from (39) reads |146j 



1 ± 




(41) 



The interaction between two electrons participating in spherically symmetric electron oscillations has been consid- 
ered in Ref. [146 . The latter predicts that there would be an effective attraction between electrons mediated by low 
energy [given by the minus sign in Eq. (41)] spherically symmetric oscillations of electrons in a quantum plasma. We 
can thus have attracting degenerate electrons in dense plasmas. The underlying physics of electron attraction here 
seems to be similar to that of the Cooper pairing of electrons in superconductors in which electrons close to the Fermi 
level attract each other due to their interactions with crystal lattice vibrations (phonon oscillations). The pairs of 
electrons act more like bosons which can condensate into the same energy level, contrary to single electrons which are 
fermions and must obey the Pauli exclusion principle. 



III. NONLINEAR EQUATIONS FOR UNMAGNETIZED QUANTUM PLASMAS 

In the preceding section, we have seen that the quantum Madelung fluid description predicts a diffraction pattern of 
a single electron or positron. However, collective interactions between an ensemble of degenerate electrons (Fermions) 
in dense plasmas are responsible for new linear and nonlinear waves and structures. 

The quantum A^-body problem is governed by the Schrodinger equation for the A^-particle wave function 
"0(91: 92, • ■ • , ^AT, t), where qj = (rj,Sj) is the coordinate (space, spin) of particle j. For identical Fermions, the 
equilibrium A^-particle wave function is given by the Slater determinant [S] 



V'(gi, 92, • ■ • , qN,t) 



1 



V'l('72,i) ^2{q2,t) 



ipNiqi^t) 
ipN(.q2,t) 



(42) 



which is anti-symmetric under odd numbers of permutations. Hence, ip vanishes if two rows are identical, which is 
an expression of the Pauli exclusion principle that two identical Fermions cannot occupy the same state. Example 
{N = 2): ip{qi,q2,t) = :^[il^i{qi,t)ip2{q2,t) - il^i{q2,t)tp2{qi,t)] so that tp{q2,qi,t) = -^plqi,q2,t) and ip{qi,qi,t) = 0. 
Due to the Pauli exclusion principle, all electrons are not permitted to occupy the lowest energy state, and in the 
ultra-cold limit when all energy states up to the Fermi energy level are occupied by electrons, there is still a quantum- 
statistical pressure determined by the Fermi pressure. 

To describe collective electrostatic oscillations in a plasma, the quantum analogue of the Vlasov-Poisson system is 
the Wigner-Poisson system, given by 



dl 
dt 



(27r)3?i4 



p2me(v— v')-A/h 



t 



A 

X A 

2' 



/(x,v',t) d^Xd-^v', 



(43) 



and 



47re 



fd v-no 



(44) 



assuming immobile ions. We note (see Appendix A) that the Wigner equation converges to the Vlasov equation for 
classical particles (electrons) when ?l — > 



dt 



V/ 



V0- 

nip ov 



(45) 



We now take the moments of the Wigner equation |43| and obtain [up to 0(h^)] the non-relativistic quantum- 
electron fluid (or the quantum Madelung fluid) equations [7J [5T] composed of the electron continuity equation 



dfie 
~dt 



V • (UeUe) = 



(46) 
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the electron momentum equation including the quantum statistical pressure and the quantum force 

me(|-+Ue-V)ue = eV(/)-— VPe + Fcj, (47) 



where (j) is determined from the Poisson equation 

V^^ = 47re(ne - Tio)- (48) 

For the degenerate Fermi-Dirac distributed plasma, one has (up to constants of order unity) the quantum statistical 
pressure for the electrons 

Pe = f ^ V""'''" , (49) 

3 V"o/ 

where D is the number of degrees of freedom in the system. It should be noted that Eliasson and Shukla [53] and 
Tsintsadze and Tsintsadze |19j have obtained different expressions for in the non-zero limit of the electron Fermi 
temperature. 

The quantum force [2(J| due to electron tunneling through the Bohm potential is 



where (f>B represents the Bohm potential. We note that the a-term [cf. (14)] does not appear in (50) due to consid- 
eration of the term up to 0{h^) in the expansion parameter, as it also happens when the mean-field approximation 
[SU] is used in deducing Fq. 

A. Nonlinear Schrodinger-Poisson equations 

By introducing the wave function 

V'(r,t) = ^ Heir, t) exp(z(pe(r, (51) 

where S is defined according to TOgUg = V^Sg and Uf. = IV'Pi it can be shown that the QHD equations [e.g. Eqs. 
(46)-(48)] are equivalent to the generahzed NLS-Poisson system [Tj [21] 



+ ^^VV + e4>^ - "^1^1^/^^ = 0, (52) 

and 



dt 2me 2nl 



V20 = 47re(|V^2| -no). (53) 

The derivation of (52) required the electron plasma flow velocity to be curl free everywhere in a singly-connected 
region, except at points where the electron number density vanishes. This is obviously not valid in general. Similar 
to quantum-fluid treatment, one should include the generalized electron vorticity D,^ = V x (|-0pUe)/|V'P, which is 



non- vanishing [see the discussions below Eq. (25)]. Equation (52) captures the two main properties of a quantum 



plasma, namely the quantum statistical pressure (fully nonlinear) and quantum dispersion effects, and is coupled self 



consistently to the electrostatic potential via the Poisson equation ( 53 1 . We thus have a nonlocal nonlinear interaction 



between the electron density and the electrostatic potential. Furthermore, we note that one-dimensional version of 



Eq. ( 52 1 without the (fhterm has also been used to describe the behavior of a Bose-Einstein condensate |147j in the 
absence of the confining potential. 

Linearization of the NLS-Poisson Equations yields the frequency of the EPOs [H [21 [H [5] 

-'^=[-pe + k'Vl + l^) ■ (54) 

One can identify two distinct dispersive effects from (54): One long wavelength regime, Vpe 3> 'hk/2me, and one short 
wavelength regime, Vpe ^ fik/2'me, separated by a critical wavenumber kcrit = ^n/Xcrit = 7r?i/meVpe 



n, 



-1/3 
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B. Inclusion of the ion dynamics 

The dynamics of low-phase speed (in comparison with the electron Fermi speed) nonlinear electrostatic ion oscilla- 
tions in a quantum electron-ion plasma is governed by the inertialess electron equation of motion 



= e\/(j) ~ — VPe 



(55) 



the ion continuity 



dt 



V • {UiVLi) = 0, 



(56) 



the ion momentum equation 



mini 



d_ 
di 



u, • V u 



-ZiTiieVcj) 



(57) 



and the Poisson equation 



= 47re(ne - ZiUi), 



(58) 



where Ui is the ion number density, is the ion fluid velocity, and Zi is the ion charge state. 

In the quasi-neutral approximation, viz. rii — ZiUi — n, we can combine Eqs. (55 1 and (57 1 to obtain 



— + • V 1 Ui = 



ZiTli 

no 



(D+2)/D 



2mpm, 



(59) 



Equations ( 56 1 and ( 59 1 are the desired set for studying nonlinear ion waves |6H 162] in a dense quantum plasma. 



IV. LOCALIZED ELECTROSTATIC EXCITATIONS 



We are now in a position to discuss nonlinear properties and the dynamics of localized electrostatic excitations in 
a dense quantum plasma based on the nonlinear equations we have developed in the preceding section. 



A. Dark solitons and vortices associated with EPOs 



In this subsection, we shall discuss the formation of one-dimensional quantized dark solitons and two-dimensional 
(2D) quantized vortices associated with the EPOs in a dense quantum plasma at nanoscales (of the order of Vpe/'^pe)- 
We note that the dynamics of dark solitons in the 2D nonlinear Schrodinger equation in a defocusing medium has been 
studied by Ivonin |148j and Ivonin et al. |149j . However, for studying the formation and dynamics of electrostatic 
nanostructures in a dense quantum plasma, we shall use the nonlinear Schrodinger-Poisson equations ^22j 



at 



and 



1*1 



(60) 



(61) 



where the normalized wave function ^' is '^j y/n^, the normalized kp is e(j)/kBTpe, and the time and space variables 



are in units of h/kBTpe and Vpe/^^pe, respectively. We have denoted A = ATrniee^ fh^riQ^^ (see Ref. [H]). 



The system of equations(|60| and (|61| is supplemented by the Maxwell equation 

iA (W^-* - «'*V*) , 



5E 

1h 



(62) 



where the electric fleld is E = 



13 



The system (60)-(62l has the following conserved integrals: the number of electrons 



N 



the electron momentum 



the electron angular momentum 



and the total energy 



v|/*r X V*c?^x, 



2+4/D 



D/{2 + D)]d^x. 



(63) 



(64) 



(65) 



(66) 



The conserved quantities are required to check the accuracy of the numerical integration of (60) and (61). 

For quasi-stationary one-dimensional structures moving with a constant speed vq, one can find localized, solitary 
wave solutions by introducing the ansatz "if = W{£,) exp{iKx — iUt), where Vl^ is a complex- valued function of the 
argument ^ = x — v^t, and K and f2 are a constant wavenumber and the frequency shift, respectively. By the choice 
K — vo/2A, the coupled system of equations can be written as 



(PW_ 



and 



where X = Q./A- 
we determine A 



ipW 
~A 



= 0, 



(67) 



(68) 



Uq/4A^ is an eigenvalue of the system. From the boundary conditions |W^| = 1 and = at |^| = co 



l/A and = 1 -I- Vq/AA. The system of Eqs. (67 1 and (68 1 admits a first integral in the form 



n = A 



dW 



d^ 



(69) 



where the boundary conditions = 1 and — Q 'At |^| = oo have be en emplo yed . 

Figure [2] shows profiles of \W\^ and (/? obtained numerically from (67l and (68I for a few values of A, where W 
was set to —1 on the left boundary and to -1-1 on the right boundary, i.e. the phase shift is 180 degrees between 
the two boundaries. The solutions are in the form of dark solitons, with a localized depletion of the electron density 
Ne = associated with a localized positive potential. Larger values of the parameter quantum coupling parameter 

A give rise to larger- amplitude and wider dark solitons. The solitons localized '"shoulders" on both sides of the density 
depletion. 

Numerical solutions of the time-dependent system of Eqs. ( 60 ) and ( 61 1 is displayed in Fig. [3] with initial conditions 
close (but not equal) to the ones in Fig. |2] Two very clear and long-lived dark solitons are visible, associated with a 
positive potential of (yS w 3, in agreement with the quasi-stationary solution of Fig. |2]for A = b. In addition there are 
oscillations and wave turbulence in the time-dependent solution presented in Fig. [3| Hence, the dark solitons seem 
to be robust structures that can withstand perturbations and turbulence during a considerable time. 

For the two-dimensional {D — 2) system, it is possible to find vortex structures of the form = tp{r) exp{is6~i^lt), 
where r and 9 are the polar coordinates defined via x = rcos{9) and y = rsm(9), is a constant frequency shift, and 
s = 0, ±1, ±2, . . . for different excited states (charge states). The index s is also known as the circulation number 
[51] . With this ansatz, Eqs. (60 1 and (61) can be written in the form 



n + A 



dr^ 



1 d 
r dr 



^ +^-1^1' 



(70) 



and 



\ dr^ r dr J ' 



(71) 
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FIG. 2: The electron density \W\'^ (the upper panel) and electrostatic potential (p (the lower panel) associated with a dark 
soliton supported by the system of equations ( |67[ ) and (68), for A = 5 (solid lines), A — 1 (dashed lines), and A = 0.2 
(dash-dotted hne). After Ref. [22]. 




FIG. 3: The time-development of the electron density l^'p (left-hand panel) and electrostatic potential (p (the right- 
hand panel), obtained from a simulation of the system of equations (601 and (611. The initial condition is 4' = 0.18-1- 
tanh[20sin(j:/10)] ex.p{iKx), with K = vo/2A, A = 5 and vq = 5. After Ref. 22_. 



respectively, where the boundary conditions ip = I and ip — dip/dr = at r = oo determine the constant frequency 
f2 = 1. Different signs of charge state s describe different rotation directions of the vortex. For s ^ 0, one must have 
"0 = at r = 0, and from symmetry considerations one has dtp/dr = at r = 0. Figure [4] shows numerical solutions 
of Eqs. ([70| and (71 1 for different s and for A = b. Here the vortex is characterized by a complete depletion of the 



electron density at the core of the vortex, and is associated with a positive electrostatic potential. 

Figure [5] shows time-dependent solutions of Eqs. (60) and (61 1 in two space dimensions for singly charged (s — ±1) 
vortices, where, in the initial condition, four vortex-like structures were placed at some distance from each other. The 
initial conditions were such that the vortices are organized in two vortex pairs, with si — -1-1, S2 = —1, S3 = —1, and 
S4 = -|-1, seen in the upper panels of Fig. [Sj The vortices in the pairs have opposite polarity on the electron fluid 
rotation, as seen in the in the upper right panel of Fig. [5j Interestingly, the "partners" in the vortex pairs attract each 
other and propagate together with a constant velocity, and in the collision and interaction of the vortex pairs (see 
the second and third pairs of panels in Fig. |5]), the vortices keep their identities and change partners, resulting into 
two new vortex pairs which propagate obliquely to the original propagation direction. On the other hand, as shown 
in Fig. |6j vortices that are multiply charged {\sj\ > 1) are unstable. Here the system of Eqs. (60 1 and (61) was again 
solved numerically with the same initial condition as the one in Fig. [5] but with doubly charged vortices si = -1-2, 
S2 = —2, S3 — —2, and S4 = -1-2. The second row of panels in Fig. [6]reveals that the vortex pairs keep their identities 
for some time, while a quasi one-dimensional density cavity is formed between the two vortex pairs. At a later stage, 
the four vortices dissolve into complicated nonlinear structures and wave turbulence. Hence, the nonlinear dynamics 
is very different between singly and multiply charged solitons, where only singly charged vortices are long-lived and 
keep their identities. This is in line with previous results on the nonlinear Schrodinger equation, where it was noted 
that vortices with higher charge states are unstable |149j . 
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FIG. 4: The electron density j^'p (upper panel) and electrostatic potential ip (lower panel) associated with a two-dimensional 
vortex supported by the system (701 and (71 1, for the charge states s = 1 (solid lines), s = 2 (dashed lines) and s = 3 
(dash-dotted lines). We used ^ = 5 in all cases. After Ref. [2^ . 



B. Localized ion wave excitations in quantum plasmas 

In his classic paper, Haas et al. developed both small and large amplitude theories for one-dimensional solitary 
ion waves in a dense quantum plasmas. They found that the dynamics of small amplitude solitary waves is governed 
by the Kortweg de-Vries (k-dV) equation 



d_U_ 



2U 



,dU_ 



1 



= 0, 



(72) 



where U represents the relative (with respect to uq) ion density perturbation, the time and space variables are in 
units of the ion plasma period cj^^/ and the electron Thomas-Fermi radius (fcBTFe/47rnoe^)^/^, respectively. 

The K-dV equation admits both the solitary and the periodic (cnoidal) waves [150] . Introducing the wave form 
U{ri = ^ — Mst), we can write (72) in the stationary frame as 



1 - 



8 ; drj^ 



- MsU + + C = 0, 



(73) 



where Mg represents the Mach number, and C is a constant of integration. In the special case when U and its 
derivatives tend to zero at ±oo, C = 0. Multiplying (73) by dU/drj one can integrate once the resultant equation, 
and express it in the form of an energy integral [1511 11521, 1153] . The resulting solitary wave solution of (73) is 



U = f/„sech^ (77/770), 



(74) 



where [/^ = (3MJ2) and tjo = {2/Ms)^/^{l - H/Sfl'^ are the maximum amplitude and the width of the soliton. We 
see that compressive solitary wave solutions are possible if < < 8. 

There also exist possibility of one-dimensional large amplitude localized ion wave excitations. To demonstrate 
this, we assume that the quantum force acting on the electrons dominates over the quantum statistical pressure, viz. 
kBTpeUe <C {h^ /4:me)d^ne/dx'^ . Hence, the electron density is obtained from [62^ 



2771, 



(75) 



The electrons are coupled with ions through the space charge electric field (— V0). 
The dynamics of singly charged ions is governed by the ion continuity 



drij 

at 



dirijUj) 
dx 



0, 



(76) 
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FIG. 5: The electron density (left panel) and an arrow plot of the electron current i (^V^** — 4'*V4') (right panel) 

associated with singly charged (|s| = 1) two-dimensional vortices, obtained from a simulation of the time-dependent system of 
equations ( |60[ ) and ( |61[ ), at times t — t = 3.3, t = 6.6 and t — 9.9 (upper to lower panels). We used A = 5. The singly 
charged vortices form pairs and keep their identities. After Ref. [22| . 



and ion momentum equation 



VTii 



d_ 
dt 



dx 



(77) 



where Ui is x component of the ion fluid velocity perturbation. The system of Eqs. (75 1-(77) is closed by the Poisson 
equation 



92 

dx 



2 = 47re(ne - n^). 



(78) 



We now look for stationary nonlinear ion wave structures moving with a constant speed uq. Hence, all unknowns 
depend only on the variable ^ = x — uqI. Defining y/n^ = V'l Eq. (75) takes the form 



+ e(pip = 0. 



(79) 



Equations (76) and (77) can be integrated once with the boundary conditions = uq and = at ^ = |c»|, and 
the results can be combined to have 



(80) 
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FIG. 6: The electron density Iv&l^ (left panel) and an arrow plot of the electron current i (5'V^'* — \['*V\&) (right panel) 
associated with double charged (|s| — 2) two-dimensional vortices, obtained from a simulation of the time-dependent system of 
Eqs. 1 60 1 and (611, at times t — 0, t — 3.3, t — 6.6 and t — 9.9 (upper to lower panels). We used A = 5. The doubly charged 



vortices dissolve into nonlinear structures and wave turbulence. After Ref. 1221 



Inserting (80) into (78 1 we obtain 



UqUo 



(81) 



Equations (79 1 and (81 1 are the desired equations for studying the nonlinear ion waves in dense quantum plasmas. 

It is convenient to introduce dimensionless quantities [see below Eq. (83l] into Eqs. (79) and (81), and rewrite 
them as 



and 



ax2 



-0, 



M 

y/M - 2$ 



= 0, 



(82) 



(83) 



where we have normalized the space variable as X = fcg^, the electron wave function as ^' = y/noijj, and the potential as 
<J> = e(j)/mic'^. Here Cq = Wpi/kg is the quantum ion wave speed and kg = {2mf,L0pf,/K)^/'^ is the quantum wavenumber. 
The quantum "Mach number" is defined as M = uo/cq. 



We note that the coupled Eqs. (82) and (83) admit a conserved quantity 



dX 



1 

2 (dX 



<^^^ - M{\/m^ - 2$ - M) = 0, 



(84) 
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where we have used the boundary conditions $ = d^/dX = d^/dX 
solitary ion wave structure, we can assume that $ = <i>,nax and ^' = 



Hence, at X = Eq. (84) yields 



0, and = 1 at \X\ 
X, as well as d'^ /dX - 



= oo. For a symmetric 
d<^/dX = at X = 0. 



^max *max 



A/) 



0. 



In the wave-breaking limit, where M = (2<i>max)^^^, we find that — <i>max^max + 2*i'max = 0, or ^'max 
the electron density will locally rise to twice the background density at wave breaking. 



(85) 

2. Accordingly, 




FIG. 7: The profiles of the potential $ (top panel) and the electron density (bottom panel) as a function of X, for different 
values of the Mach number: M = 1.1 (dashed curves), M = 0.9 (solid curve), and M — 0.75 (dotted curve). After Ref. |62| . 



In Fig. [t] numerical solutions of Eqs. (82 1 and (83) are displayed, showing profiles of the electrostatic potential and 
electron number densities for different values of M. We see that both the electrostatic potential and electron density 
have localized and strongly peaked maxima and an oscillatory tail. The latter is in sharp contrast to the classical 
(non-quantum) case, where the ion acoustic solitary waves have a monotonic profile, which in the small amplitude 
limit, where the system is governed by the Korteweg-de Vries equation, assumes a secant hyperbolicus shape. We 
observe from Fig. |7]that M = 1.1 is close to the wave breaking limit above which there do not exist solitary wave 
solutions. Our numerical investigation also suggests that there is a lower limit of M (slightly lower than 0.75), below 
which the solitary wave solution vanishes. 



V. QUANTUM FLUID TURBULENCE 



The statistical properties of turbulence and its associated electron transport at nanoscales in quantum plasmas 
has been investigated in both 2D and 3D by means of the coupled NLS and Poisson equations [66j |67]. It has 
been found that the nonlinear coupling between the EPOs of different scale sizes gives rise to small-scale electron 
density structures, while the electrostatic potential cascades towards large-scales. The total energy associated with 
the quantum electron plasma turbulence, nonetheless, processes a characteristic, non-Kolmogorov-like spectrum. The 
electron diffusion caused by the electron fluid turbulence is consistent with the dynamical evolution of turbulent mode 
structures. 

To investigate the quantum electron fluid turbulence in 3D, we use the nonlinear Schrodinger- Poisson equations 

EH El [67] 



HV'^^ + if^ - ^ 0, 



and 



1, 



(86) 



(87) 



which govern the dynamics of nonlinearly interacting EPOs of different wavelengths. In Eqs. (86) and (87) the wave 



function is normalized by y^^o, the electrostatic potential by ksTp^/e, the time t by the electron plasma period uj 



and the space r by the Thomas-Fermi Debye radius Vpe/i^pe- We have introduced the notation ^/H = hujpe/V^kBTp, 
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The nonlinear mode coupling interaction studies are performed to investigate the multi-scale evolution of a decaying 
3D electron fluid turbulence, which is described by Eqs. ([SG]) and (87 1. All the fluctuations are initialized isotropically 



(no mean fields are assumed) with random phases and amplitudes in Fourier space, and evolved further by the 
integration of Eqs. (86 1 and (87), using a fully de- aliased pseudospectral numerical scheme |154j based on the Fourier 
spectral methods. The spatial discretization in our 3D simulations uses a discrete Fourier representation of turbulent 
fluctuations. The numerical algorithm employed here conserves energy in terms of the dynamical fluid variables 
and not due to a separate energy equation written in a conservative form. The evolution variables use periodic 
boundary conditions. The initial isotropic turbulent spectrum was chosen close to fc~^, with random phases in all 
three directions. The choice of such (or even a flatter than —2) spectrum treats the turbulent fluctuations on an equal 
footing and avoids any influence on the dynamical evolution that may be due to the initial spectral non-symmetry. 
The equations are advanced in time using a second-order predictor-corrector scheme. The code is made stable by a 
proper de-aliasing of spurious Fourier modes, and by choosing a relatively small time step in the simulations. Our 
code is massively parallelized using Message Passing Interface (MPI) libraries to facilitate higher resolution in a 3D 
computational box, with a resolution of 128'^ grid points. 

We study the properties of 3D fluid turbulence, composed of nonlinearly interacting EPOs, for two specific physical 
systems. These are the dense plasmas in the next generation laser-based plasma compression (LBPC) schemes |;81j as 
well as in superdense astrophysical objects [96l|97l[98] (e.g. white dwarfs). It is expected that in LBPC schemes, the 
electron number density may reach 10^^ cm~'^ and beyond. Hence, we have ujr>e = 1.76 x 10^^ s-\ Tp = 1.7 xlO'^ 
hujpe = 1.7 X 10~^ erg, and H = 1. The Fermi Debye length Xd 



■pe — ^.,^^^Kj a ,j.F = i.^xiu - erg, 
0.1 A. On the other hand, in the interior of white 



yielding LUp^ = 5.64 x lO^'' s'^, Tp = 1.7 x 10"^ erg 



dwarfs, we typically have hq 

H « 0.3, and Xd = 0.025 A. The numerical solutions of Eqs. ([sej) and ([87| for H_ 
to no 



— 10^'' cm and uq ~ 



Tiujpe — 5.64 X 10 ^ erg, 
0.025 (corresponding 



- ^ 1 and 7J 

respectively) are displayed in Figs, [s] and [9] respectively, which are the 
electron number density and electrostatic (ES) potential distributions in the (x, 2/)-plane. 



(a) \vf, H=0.4 (b) * 




FIG. 8: Small scale fluctuations in the electron density resulted from a steady turbulence simulations of our 3D electron plasma, 
for H — 0.4 (top panels) and H = 0.01 (bottom panels). Forward cascades are responsible for the generation of small-scale 
fluctuations seen in panels (a) and (c). Large scale structures are present in the electrostatic potential, seen in panels (b) and 
(d), essentially resulting from an inverse cascade. After Ref. [67) . 

Figures [8] and [9] reveal that the electron density distribution has a tendency to generate smaller length-scale struc- 
tures, while the ES potential cascades towards larger scales. The co-existence of the small and larger scale structures 
in turbulence is a ubiquitous feature of various 3D turbulence systems. For example, in 3D hydrodynamic turbulence, 
the incompressible fiuid admits two invariants, namely the energy and the mean squared vorticity. The two invari- 
ants, under the action of an external forcing, cascade simultaneously in turbulence, thereby leading to a dual cascade 
phenomena. In these processes, the energy cascades towards longer length-scales, while the fiuid vorticity transfers 
spectral power towards shorter length-scales. Usually, a dual cascade is observed in a driven turbulence simulation, in 
which certain modes are excited externally through random turbulent forces in spectral space. The randomly excited 
Fourier modes transfer the spectral energy by conserving the constants of motion in fc-space. On the other hand, in 
freely decaying turbulence, the energy contained in the large-scale eddies is transferred to the smaller scales, leading 
to a statistically stationary inertial regime associated with the forward cascades of one of the invariants. Decaying 
turbulence often leads to the formation of coherent structures as turbulence relaxes, thus making the nonlinear in- 
teractions rather inefficient when they are saturated. The power spectrum exhibits an interesting feature in our 3D 
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FIG. 9: Power spectrum of 3D EPOs in the forward cascade regime. A Kolmogorov-like spectrum ~ k ^^^^ is observed for 
H = 0.4. The spectral index changes as a function of H. The numerical resolution is 128^. After Ref. [67j . 



electron plasma system, unlike the 3D hydrodynamic turbulence |155l 11561 11571 [T55] . The spectral slope in the 3D 

fc-3/2, 



quantum electron fluid turbulence is close to the Iroshnikov-Kraichnan power law |159l 1160] k^'^' , rather than the 
usual Kolomogrov power law |155j k~^^^. We further find that this scaling is not universal and is determined critically 
by the quantum tunneling effect. For instance, for a higher value of H— 1.0 the spectrum becomes more flat (see Fig|9|. 
Physically, the flatness (or deviation from the k~^/^), results from the short wavelength part of the EPOs spectrum 
which is controlled by the quantum tunneling effect associated with the Bohm potential. The peak in the energy 
spectrum can be attributed to the higher turbulent power residing in the EPO potential, which eventually leads to 
the generation of larger scale structures, as the total energy encompasses both the electrostatic potential and electron 
density components. In our dual cascade process, there is a delicate competition between the EPO dispersions caused 
by the statistical pressure law (giving the k^Vp term, which dominates at longer scales) and the quantum Bohm 
potential (giving the h'^k'^/4ml term, which dominates at shorter scales with respect to a source). Furthermore, it is 
interesting to note that exponents other than k~^^^ have also been observed in numerical simulations [16111162] of the 
Charney and 3D incompressible Navier-Stokes equations. Recently, Paoletti et al |163l have examined the velocity 
statistics of quantum turbulence in supcrfluid "^He, and found that it significantly differs from classical turbulence 
due to the topological interactions of vortices that are different from those in classical fluids. 
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FIG. 10: Time evolution of an effective electron diffusion coefficient associated with the large-scale electrostatic potential and 
the small-scale electron density, for H = 0.4, H = 0.1 and H = 0.01. Smaller values of H corresponds to a small effective 
diffusion coefficient, which characterizes the presence of small-scale turbulent eddies that suppress the electron transport. After 
Ref. [ST]. 



21 



We finally estimate the electron diffusion coefficient in the presence of small and large scale turbulent EPOs in our 
quantum plasma. An effective electron diffusion coefficient caused by the momentum transfer can be calculated from 
£»e// = /o^(P(r,i) •P(r,i + i')>'^^'> where P is electron momentum and the angular bracket denotes spatial averages 
and the ensemble averages are normalized to unit mass. The effective electron diffusion coefficient, D^ff, essentially 
relates the diffusion processes associated with random translational motions of the electrons in nonlinear plasmonic 
fields. We compute D^ff in our simulations, to measure the turbulent electron transport that is associated with the 
turbulent structures that we have reported herein. It is observed that the effective electron diffusion is lower when 
the field perturbations are Gaussian. On the other hand, the electron diffusion increases rapidly with the eventual 



formation of longer length-scale structures, as shown in Fig. 10 The electron diffusion due to large scale potential 
distributions in quantum plasmas dominates substantially, as depicted by the solid-curve in Fig. |10[ Furthermore, 
in the steady-state, nonlinearly coupled EPOs form stationary structures, and i^e// saturates eventually. Thus, 
remarkably an enhanced electron diffusion results primarily due to the emergence of large-scale potential structures 
in our 3D quantum plasma. 

VI. KINETIC PHASE-SPACE STRUCTURES 

In the preceding sections, we have discussed the properties of quantized coherent structures and 3D quantum 
electron ffuid turbulence based on the coupled Schrodinger and Poisson equations. Thus, it has been assumed that 
nonlinearly interacting plasma waves are spontaneously created by some known physical processes (e.g. the beam- 
plasma instability) in quantum plasmas. 

The formation of electrostatic kinetic phase space structures, based on the Vlasov-Poisson equations, in classical 
plasmas has been well documented |164j . In the following, we shall discuss quasilinear aspects |165j of the EPOs that 
are governed by the Wigner-Poisson system (i.e. a quantum analogue of the Vlasov-Poisson system). Specifically, 
we focus on kinetic phase-space nonlinear structures arising from the trapping of electrons in the finite amplitude 
wave potential and the self-consistent modification of the electron distribution function in the presence of nano-kinetic 
structures. 

To study the differences in the nonlinear evolution of the Wigncr and Vlasov equations, we have simulated the well- 
known bump-on-tail instability |165j . whereby a high- velocity beam is used to destabilize a Maxwellian equilibrium. 
We use the initial condition / = (1 -|- S){no/^/2^TVth)[0■&GXJp{—v^/2vff^) + 0.4exp(— 2(w — 2.5vth)'^ /v^,J], where 6 
represents random fluctuations of order 10~^ that help seed the instability (see Fig. 1). Here vth = \/krBT^Jml is the 
electron thermal speed. We use periodic boundary conditions with spatial period L — 407rA£ie, where Xoe = Vthl^p is 
the Debye length. Three simulations were performed, with different values of the normalized Planck constant, defined 
as H — hujpe / mv^i^: H = (Vlasov), H = 1, and H = 2. 

In order to highlight the transient oscillations in velocity space, we first perturb the above equilibrium with a 
monochromatic wave having kXjjf, = 0.25 (i.e., a wavelength of STrA^e). Figure [TT] shows the results from simulations of 
the Wigner-Poisson and Vlasov-Poisson systems. In both simulations, due to the bump-on-tail instability, electrostatic 
waves develop nonlinearly and create periodic trapped-particle islands (electron holes) with the wavenumber k = 
0.25\J)],. The theory predicts the formation of velocity-space oscillations in the Wigner evolution, which should be 
absent in the classical (Vlasov) simulations. This is the case in the results presented in Fig. 2, where the oscillations 
are clearly visible. 

When the initial excitation is broad-band (i.e., wavcnumbcrs 0.05 < fcAue < 0.5 are excited), the electron holes 



start merging together at later times due to the sideband instability |166|, 1167] (see Fig. 12 1. At this stage, mode 
coupling becomes important and quasilinear theory is not capable of describing these effects. As the system evolves 
toward larger spatial wavelength, the evolution becomes progressively more classical, with the appearance of a plateau 
in the resonant region. Nevertheless, at ujpet = 500 the Wigner solution still displays some oscillatory behavior in 
velocity space, which is absent in the Vlasov evolution. 

From the experimental viewpoint, recent collective x-ray scattering observations in warm dense matter |118j revealed 
a measurable shift in the plasmon frequency due to quantum effects. In the nonlinear regime (strong excitations), this 
effect could lead to trapping of electrons in the wave potential of the plasmons, and the subsequent formation of the 
kind of phase space structures discussed here. Furthermore, we note that there also exists a theoretical description 
[168] of quantum corrected electron holes based on the perturbative treatment of the Wigner-Poisson equations. 

VII. MAGNETIC FIELDS IN QUANTUM PLASMAS 



There are several mechanisms by which magnetic fields in classical plasmas can be generated. They include i) 
non-parallel density and temperature gradients (the so called Biermann battery |169] ). ii) the electron temperature 
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FIG. 11: (Color online) Simulations of the Wigner-Poisson and Vlasov-Poisson systems, at time Up^t — 200, for a.) H — 
(Vlasov), h) H — 1 (Wigner), and c) H — 2 (Wigner). Initially monochromatic spectrum. Top panels: electron distribution 
function f{x, v) in phase space. Bottom panels: spatially averaged electron distribution function F{v) in velocity space. After 
Ref. [1651. 



anisotropy (known as the Weibel instability |170j '). ill) counterstreaming electron beams |171l I172j . and iv) the 
ponderomotive forces of laser beams [1731 [Ml UlSl [Ml UlTl HlHl [179] . 

In the following, we discuss two possibilities for the magnetic field generation in quantum plasmas. 



A. The quantum Weibel instability 



We first discuss linear and nonlinear aspects [1801 11811 1182LIT55] of the Weibel instability that is driven by equilibrium 
Fermi-Dirac electron temperature anisotropic distribution function in a nonrelativistic dense quantum plasma. It is 
well known |184| that a dense quantum plasma with an isotropic equilibrium distribution function does not admit 
any purely growing linear modes. This can be verified, for instance, from the expression for the imaginary part of 
the transverse dielectric function, as derived by Lindhard |185j . for a fully degenerate non-relativistic Fermi plasma. 
It can be proven (see Eq. (30) of |186j ) that the only exception would be for extremely small wavelengths, so that 
k > 2kp, where kp the characteristic Fermi wave number of the system. However, in this situation the wave would be 
super- luminal. On the other hand, in a classical Vlasov-Maxwell plasma containing anisotropic electron distribution 



function, we have a purely growing Weibel instability |170j . via which dc magnetic fields are created. The electron 
temperature anisotropy may arises due to the heating of the plasma by laser beams |187j . where there is a signature 
of the Weibel instability as well. In the next generation intense laser-solid density plasma experiments, it is likely that 
electrons would be degenerate and that electron temperature anisotropy may develop due to an anisotropic electron 
heating by intense laser beams via resonant absorption, similar to the classical laser-plasma-interaction case |188j . 

Consider linear transverse waves in a dense quantum plasma composed of the electrons and immobile ions, with 
k • E = 0, where k is the wave vector and E is the wave electric field. Following the standard procedure, one then 
obtains the general dispersion relation [TJ 11821 1183] for the transverse waves of the Wigner-Maxwell system 



; 2 2 

k c 



0, 



where v = {v^^Vy^Vz) is the velocity vector, and fo{vx,Vy,Vz) is the equilibrium Wigner function associated to Fermi 
systems. For spin 1/2 particles, the equilibrium pseudo distribution function is in the form of a Fermi-Dirac function. 
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FIG. 12: (Color online) Simulations of the Wigner-Poisson and Vlasov-Poisson systems, for Upet — 500, a) for H = (Vlasov) 
and h) H = 1 (Wigner). Initially broad wavenumber spectrum. Top panels: electron distribution function f{x,v) in phase 
space. Bottom panels: spatially averaged electron distribution function F{v) in velocity space. After Ref. [165] . 



Here we allow for velocity anisotropy and express 

/o = 

exp 



2 y kbT± kbT^ 

where /i is the chemical potential, and the normalization constant is 

no 



/3m 



Li3/2(-e^^) V 2tt ) ^ \2Tin) 



(89) 



(90) 



Here 'Liz/2 is a polylogarithm function |189lll90j . Also, (3 — I /[kb{T1T\\Y''^], where Tl and T|| are related to velocity 
dispersion in the direction perpendicular and parallel to z axis, respectively. In the special case when Tl — Ty, the 
usu al Fermi-Dirac equilibrium is recovered. The chemical potential is obtained by solving the normalization condition 
(90 1, yielding, in particular, ^ = Ep in the limit of zero temperature, where Ep = (37r^no)^^'^?i^/(2r7ie) is the Fermi 



energy. Also, the Fermi-Dirac distribution /(k ), where k is the appropriated wave vector in momentum space, is 
related to the equilibrium Wigner function (89 1 by /(k) — (l/2)(27r?i/TO)'^/o(v), with the factor 2 coming from spin 
[191} 1192] . However, these previous works refer to the cases where there is no temperature anisotropy. Notice that 
it has been suggested [193' that in laser plasmas the Weibel instability is responsible for further increase of T|| with 
time. 



Inserting (89 1 into (88 1 and integrating over the perpendicular velocity components, we obtain 



Til 



where 
Wq 



2^7rHUzf2{-ef'f^) J 



dv 



Lio 



|-exp -('^+ y) j-Liaj-exp -(^ 




(91) 



(92) 



In (92 1, Li2 is the dilogarithm function |189l 1190] . H = 7lfc/(meWj[) is a characteristic parameter representing the 
quantum diffraction effect, ^ = w/(fc?;||), and v = Vyjv^^. with = (2KBT\\l'ni^^l'^ . In the simultaneous limit of a 



24 



small quantum diffraction effect {H <C 1) and a dilute system {e^^^ <C 1), it can be shown that Wq 



-i-f^(e), 

where Z is the standard plasma dispersion function |194j . It is important to notice that either (88 1 or (91 1 reproduces 
the transverse dielectric function calculated from the random phase approximation for a fully degenerate quantum 
plasma |185j . in the case of an isotropic system. The simple way to verify this equivalence is to put T± = Ty in (88 1 
and then take the limit of zero temperature, so that /o 



3no/(4^y|J for 



v| < Vpe, and /o = otherwise. 
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FIG. 13: The growth rate for the Weibel instability of a dense Fermionic plasma with no = 10 m~ {u)pe = i.8 x iU'" s 
and = 5, relevant for the next generation inertially compressed material in intense laser-sohd density plasma interaction 
experiments. The temperature anisotropics are Tx/Tn — 3 (dashed line), Tl/T|| = 2 (solid hne) and T^/T\\ — 1.5 (dotted 
line), yielding, respectively, Tn = 3.9 x 10^ K, Tn = 5.2 x lO'^K and T|| = 6.3 x lO'^K. After Ref. [US]. 
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FIG. 14: The growth rate for the Weibel instability of a dense Fermionic plasma with no — 10'^"' m (uipe — 1.8 x 10^* s ^). 
Here the temperature anisotropy is Ti/Tn — 2. We used — 1 (dashed line), Pfj, = 5 (solid line) and = 10 (dotted line), 
yielding Tn = 1.6 x lO'' K, Tn = 5.2 x lO'' K and Tn = 2.6 x 10® K, respectively After Ref. 



We next solve our new dispersion relation (91 1 for a set of parameters that are representative of the next gener- 
ation laser-solid density plasma interaction experiments. The normalization condition (90 1 can also be written as 



Li3/2[-exp(/3^)] = (4/3V^)(^^;F)^/^ which is formally the same relation holding for isotropic Fermi-Dirac equi 
libria [8]. For a given value on the product f3fi and the density, this relation yields the value /3, from which the 
temperatures T± and T|| can be calculated, if we know Tj_/T\\. Consider only purely growing modes. From the 
definition (92), one can show that Wq when = 17 — > for a finite wavenumber k. From (91 1 we then obtain 

the maximum wavenumber for instability as /cmax — {^pe/c)^jT\jT\\--l. When Tl/T|| ~> 1, the range of unstable 
wavenumbers shrinks to zero. In Figs. 13 and 14 we have used the electron number density no = lO'^'^m"'^, which 



can be obtained in laser-driven compression schemes. The growth rate for different values on T^/T\\ is displayed in 

{ujpelc)yfT^/T\\ 



Fig 



13 



1, as predicted, and that the 



We see that the maximum unstable wavenumber is /cmax 

m growth rate occurs at fc « A:„iax/2. Figure 13 also reveals that the maximum growth rate of the instability 
is almost linearly proportional to T^/T\\ — 1. In Fig. 14 we have varied the product which is a measure of the 



degeneracy of the quantum plasma. We see that for p/i larger than 5, the instability reaches a limiting value, which 
is independent of the temperature, while thermal effects start to play an important role for of the order unity. 
From several numerical solutions of the linear dispersion relation, it was deduced |183j an approximate scaling law 

for the instability as ^maxl^pe = constant x nlj^ {T±^/T\\ — 1), where the constant is approximately 8.5 x 10~^'^s~^m. 
Using that no = {2meEF /fi^ f'^ / {^n^) « 1.67 x W^{Ef /meC^f/^ , we have 



= 0.10 



Wo 



Ef 



1/2 



Til 



- 1 



(93) 



for the maximum growth rate of the Weibel instability in a degenerate Fermi plasma. This scaling law, where the 
growth rate depends on the Fermi energy and the temperature anisotropy, should be compared to that of a classical 
plasma |188l I195j , where the growth rate depends on the thermal energy and the temperature anisotropy. 



25 



5000 



„ 3000 
^ 2000 



3 10 20 30 40 5( 



10 20 30 40 50 

X!L> /C 



-^-0.01 



10.005 



FIG. 15: The magnetic field components By (top panel) and (bottom panel) as a function of space and time, for (3^ = 5 
and T^/T\\ — 2. The magnetic field has been normalized by ujpeme/e. We see a nonlinear saturation of the magnetic field 
components at an amplitude of ~ 0.01. After Ref. |183] . 
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FIG. 16: The maximum of the magnetic field amplitude, B = {By + B^)^^^, over the simulation box (top panel), and 
the logarithm of the magnetic field maximum (bottom panel) as a function of time, for Tx/Tji — 2 and = 5. The 
magnetic field has been normalized by uipsnris/e. From the logarithmic slope of the magnetic field in the linear regime one finds 
7 « Aln(Bmax)/At « Q.Oltjp. After Ref. [15^ . 



For a Maxwellian plasma, it has been found |196| that the Weibel instability saturates nonlinearly once the magnetic 
bounce frequency ujc = eB/m^c has increased to a value comparable to the linear growth rate. In order to assess 
the nonlinear behavior of the Weibel instability for a degenerate plasma, we have carried out a kinetic simulation of 
the Wigner-Maxwell system. We have assumed that the quantum diffraction effect is small, so that the simulation 
of the Wigner equation can be approximated by simulations of the Vlasov equation by means of an electromagnetic 
Vlasov code |197j . As an initial condition for the simulation, we used the distribution function (89 1. In order to 
give a seed for any instability, the plasma density was perturbed with low-frequency fluctuations (random numbers) . 
The results are displayed in Figs. [l5|and 16 for the parameters /3/x = 5 and T±/T\\ — 2, corresponding to the solid 

^"TshowsTne magnetic field components as a function of space and time. We see 



lis and 16 for the parameters /3/x 
lines in Figs. [13] and [14] Figure [^ 
that the magnetic field initially grows, and saturates to steady state magnetic field fiuctuations with an amplitude of 
eB /mecujpe w 0.008. The maximum amplitude of the magnetic field over the simulation box as a function of time is 
shown in Fig. 16 where we see that the magnetic field saturates at eB/nieUjpe ~ 0.0082, while the linear growth rate 



of the most unstable mode is jmax/^pe ~ 0.009. Similar to the classical Maxwellian plasma case |196j . we can thus 
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estimate the generated magnetic field as 



^^meCTW^ (94) 



for a degenerate Fermi plasma. For our parameters parameters relevant for intense laser-solid interaction experiments, 
we will thus have magnetic fields of the order 10^ Tesla (one gigagauss). 

B. Dense plasma magnetization by the electromagnetic wave 

The second example of the dense plasma magnetization can occur in the presence of streaming electrons and large 
amplitude electromagnetic waves. Here, the nonstationary ponderomotive force of the electromagnetic wave would 
create slowly varying electric fields and currents, which generate d.c. magnetic fields in dense plasmas. 

Let us consider the propagation of the electromagnetic wave with the electric field E(r, t) = (l/2)Eo(a;, t) exp{—iujt+ 
ikx) + c.c, in an unmagnetized non-relativistic dense plasma with streaming electrons (with the drift velocity uz, 
where u is the magnitude of the electron drift speed and z is the unit vector along the z axis in a Cartesian coordinate 
system; typically u is much smaller than the electron Fermi speed) and immobile ions. Here, Eo(a;,i) is the envelope 
of the electromagnetic field at the position r and time t, and c.c. stands for complex conjugate. The frequency uj and 
the wave vector k = fcic, where x is the unit vector along the x axis, are related by |198j 

where N is the index of refraction and fig — Tik'^ /2m^, 

The electromagnetic wave exerts a ponderomotive force Fp = Fps+Fpt on the plasma electrons, where the stationary 
and non-stationary ponderomotive forces |173j are, respectively, 

Fps = ^^^^V|Eop, (96) 

and 

1 k d[u;HN~l)]d\Eo\' 
The ponderomotive force pushes the electrons locally, and creates the slowly varying electric field 



1 rIA 1 

E« = -V(/),--^ = — Fp, (98) 

C Ot TIqC 



where the scalar and vector potentials are, respectively, 

'"-''lEol", (m) 



and 

c k d[uj^{N-l)]. 



The induced slowly varying magnetic field is then = V x A^. Noting that 



For (100) 



d[u,\N~l)] ^ 2^fcVo.^, 
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we can express the magnitude of the magnetic field as 



IB, 



(102) 



where L is scale length of the envelope |Eop and — {k'^Vp^ 



f]2)i/2 ^ (7ifcV2me)[l+4(37r2no)2/VA:2]i/2^ We note 

from (102) that the magnetic field strength is proportional to u^, and attains a large value when a; ~ fi. The electron 

gyrofrequency flc is 



elBJ k^V^u^uj 



^c- — = ,U 2 (103) 



where Vo = e|Eo|/mea; is the electron quiver velocity in the electromagnetic field. 



VIII. DYNAMICS OF ELECTROMAGNETIC WAVES IN DENSE PLASMAS 



We here consider various electromagnetic (EM) wave modes and their nonlinear interaction in dense plasmas. As 
examples, we will discuss spin waves, which occur in a magnetized plasma due to the spin-1/2 effect of electrons and 
positrons. We also consider nonlinear interactions between finite amplitude electromagnetic and electrostatic waves 
in dense plasmas. We focus on the underlying physics of stimulated scattering instabilities and the wave localization 
due to the parametric interactions involving the radiation pressure |178] . 



A. Electromagnetic spin waves in magnetized plasmas 

In the presence of an external magnetic field, the quantum description of the linear kinetics of a dense coUisionless 
plasma is complicated due to quantization of the gyromagnetic motion and the inclusion of the electron and positron 
spin-1/2. In the past, many authors |199l 12001 12011 [2U^ I203j investigated the high-frequency conductivity, longitudinal 
and transverse dielectric responses in dense magnetized plasmas. The propagation characteristics of high-frequency 
electromagnetic waves in quantum magnetoplasmas is different from that in a classical magnetoactive plasma |204j . 

Recently, Oraevsky et al. [7(1 found a new electromagnetic spin wave whose electric field is parallel to zBq, and 
which propagates across z, where z is the unit vector along the z axis in a Cartesian coordinate system and is the 
strength of the ambient magnetic field. The spin wave accompany the magnetization current due to spinning electron 
motion and the wave frequency is obtained from |70] 



k 



(104) 



where k± is the perpendicular component of the wave vector k — xkj_, x is the unit vector transverse to z, v,, is the 
electron collision frequency, nd ^ rrieC^. The third term in the right-hand side of (102) represents the electron 
spin-magnetic resonance at the electron gyrofrequency uoce = eBQ/m^c. 
Assuming that /c^c^ -I- ^ w^g, we obtain from (102) 



27r^|no 



klc^ 



klc^+UJ^^{l-iu,/Lu)~UJ^ 



(105) 



For <C w 



the damping rate of the spin mode is 



Im'^ - , ^ „2.'2 Tu2,2 , ..2 ^rT2- (106) 



The ponderomotive force of the spin wave can create compressional magnetic field perturbation due to an inverse 
Cotton-Mouton/Faraday effect [205^ . 
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B. Nonlinear ly coupled EM waves 

Finite amplitude electromagnetic waves in quantum magnetoplasmas interact nonlinearly among themselves. In 
this subsection, we shall use the generalized Q-MHD equations to obtain compact nonlinear equations for the electron 
magnetohydrodynamic (EMHD) and the Hall-MHD plasmas, and show how does the density, the fluid velocity, and 
magnetic field perturbations are coupled in a non-trivial manner. 

The governing nonlinear equations for the electromagnetic waves in dense magnetoplasmas are the quantum mag- 
netohydrodynamic equations composed of the continuity equation 

— ^ + V • K,.Ue,.) = 0, (107) 

at 

the electron and ion momentum equations, respectively, 

neirie (^^+Ue- V j Ue = -n^e + ^Ue X - VPe + Fge, (108) 

niTTii Ui • Ui = ZiCUi + x , (109) 

the Faraday law 

cVxE = — — , (110) 
the Maxwell equation including the magnetization spin current 

An 1 aE 

VxB= — (Jp-|-J„0 + - — , (111) 

C COT 

where the pressure for non-relativistic degenerate electrons reads |206] 

3/2' 



3(27r)2?i2c 



(112) 



where 77.^ = 0, 1, 2, rimaxi and the value of n„iax is fixed by the largest integer nj^ that satisfies Ep — nj^tkUi-e ^ 0. 
The sum of the quantum Bohm and intrinsic angular momentum spin forces is 



In Eqs. ( 107 )-( 1131, rij is the number density of the particle species j (j equals e for the electrons, and i for the ions) , 
Uj is the particle fluid velocity, and B — |B|. We have denoted the plasma current density 3p — —n^eu,, + ZiUieUi 
and the electron magnetization spin current density 3m = V x M, where the magnetization for dynamics on a time 
scale much slower than the spin precession frequency for an electron Fermi gas reads pUBj M — {ne^i%/kBTpe)'B. 

1. Nonlinear EMHD 

First, we present the generalized nonlinear EMHD equations for a dense magnetoplasma. Here the ions form the 
neutralizing background. The wave phenomena in the EMHD plasma will occur on a time scale much shorter than 
the ion plasma and ion gyroperiods. In equilibrium, we have 

neo = ZiTHo = no. (114) 
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The relevant nonlinear EMHD equations are 



+ V-Kue)=0, (115) 



the electron momentum equation (1081, Faraday's law ( |110 l, and the electron fluid velocity given by 



^ Jm c(V X B) ^ 1 

eUe Airerie eUf. dt 

We observe that the quantum tunneling and spin forces play an important role if there are slight electron density and 
magnetic field inhomogeneities in dense plasmas. The nonlinear EMHD equations are useful for studying collective 
electron dynamics in metallic and semiconductor nanostructures |112j . 

2. Nonlinear Hall-MHD 

Second, we derive the modified nonlinear Hall-MHD equations in a dense electron-ion plasma. The Hall-MHD 
equations shall deal with the wave phenomena on a time scale larger than the electron gyroperiod. The relevant 
nonlinear Hall-MHD equations are the electron and ion continuity equations, the inertialess electron momentum 
equation 

E-K-UexB-^-?^ = 0, (117) 



Faraday's law (110), the ion momentum equation 



nt-m., ( ^ Uj • V j = Z,en, -|- ^u^ x B ) , (118) 



and the electron fluid velocity given by 



_ c(V X B) 



(119) 



where we have neglected the displacement current, since the Hall-MHD plasma deals with electromagnetic waves 
whose phase velocity is much smaller than the speed of light in vacuum. 



We now eliminate the electric field from (1181 by using (117), obtaining 



Tiimi 



dt 



u,: • V u,: = Zierii 



-(Ui - Ue) X B 

c riee ripe 



Furthermore, eliminating Up from (120) by using (119), we have 



(120) 



mrrii 



u, 



= ZiCrii i - 



c(V X B) 



in 



with the quasi-neutrality condition Up = ZiUi 



Finally, by using (117) we can eliminate E from (108), obtaining 



X B 



VPp 



Hpe 



(121) 



dB 

~dt 



V X 



(u, X B) . 



Jm X B 

eZiU., 



47r 



(V X B) X B 



(122) 



The ion continuity equation, Eqs. (121 1 and (122), together with (118) and ZiUn — Upi, where Upi^n <^ neo,io are 
the desired generalized nonlinear equations for the low-frequency (in comparison with the electron gyrofrequency) , 
low-phase velocity (in comparison with the speed of light in vacuum) in a Hall-MHD dense plasma. They can be 
used to investigate the multi-dimensional linear and nonlinear waves (e.g. magnetosonic solitons |209p . as well as 
nanostructures and turbulences [210j in dense quantum magnetoplasmas. 
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C. Stimulated scattering instabilities 



Nonlinear interactions between the high-frequency EM waves and low-frequency electrostatic waves give rise to 
stimulated scattering instabilities in classical plasmas |130|, I21H \212\ 1213] . There also exists possibility of exciting 
plasma waves by the high-frequency EM waves in quantum plasmas due to the parametric instabilities |214j . The 
governing equations for the high-frequency electromagnetic waves |64j and the radiation pressure driven modified 
Langmuir and ion-acoustic oscillations in an unmagnetized quantum plasma are, respectively, 

— -c''V^ + ujIAa + ujI,J^^A^O, (123) 
£ + < - ^^leV^ + ^^') ^ - -^WW, (124) 

and 



'^'yky^ + j^^')--.;Ar^^'\M^ (125) 



where A is the vector potential of the high-frequency electromagnetic wave, and ni is the electron density perturbation 
of the low-frequency oscillations (viz. the modified EPOs and ion waves). 



Combining Eqs. (123l-(125l we thus simply obtain the nonlinear dispersion relations 



2 n2 _ ^S'^pe^ |Ao 

— itr, — — 



and 



^2_^2^_^I^^^^ 1 (127) 

where nR = (w^^ + 3k^V^j5 + n^k^/Amlf/^ and VIb = (fc^Cl, + fi^k^ /memif^ . 

Equations (126) and (127) admit stimulated Raman, stimulated Brillouin and modulational instabilities of the 
electromagnetic pump (with the amplitude Aq and the frequency t^o = (^qC^ + ^pe)^^^) ™ dense quantum plasmas. 
We have denoted D± — ±2a;o(ij-' — k • Vg =F 5), where Vg — c^ko/wo is the group velocity of the EM pump wave 
frequency, and 5 — k'^c^ /2ijJq is the small frequency shift arising from the nonlinear interaction of the pump with 
the electrostatic perturbations (w,k) in our quantum plasma. For stimulated Raman and Brillouin scatterings, the 
EM pump wave is scattered off the resonant electron plasma wave and resonant ion wave, respectively, while for the 
modulational instability the electron and ion plasma oscillations are off resonant, so are the EM sidebands. 

For three- wave decay interactions, we suppose that Z?_ = and ^ 0. Thus, we ignore Z?-|_ from (126) and 
(127). Letting ui = + ijn and to = flB + i^B in (126) and (127), respectively, and assuming w — k • + 5 = i^r^b, 
we obtain the the growth rates for stimulated Raman and Brillouin scattering instabilities (denoted by the subscripts 
R and B, respectively, 

2\/2u)i^^RmeC 

and 

. ek\A,? ^^^^^ 
2^/ ZLUQiiBTnemiC 

where |k-Vg — (5| ^ i^R,ilB- We note that the growth rates for stimulated Raman and Brillouin scattering instabilities 
are inversely proportional to the square roots of ^Ir and ^Ib, which depend on the quantum parameters. Stimulated 
Raman and Brillouin scattering instabilities of the EM wave off the ES waves should provide invaluable information 
regarding the density fluctuations and equation of states that might exist in dense quantum plasmas. 

The quantum corrected 3D Zakharov equations |127j have been derived by Haas |215j and Haas and Shukla |216j . 
who demonstrated that the dispersive effects associated with quantum corrections can prevent the collapse of localized 
Langmuir envelope electric fields in both two and three spatial dimensions. 
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D. Self-trapped EM waves in a quantum hole 



Nonlinear interactions between large amplitude electromagnetic waves and electrostatic plasma waves can produce 
nonlinear nanostructures composed of a density cavity that traps the electromagnetic wave envelope. Here we demon- 
strate trapping of intense electromagnetic waves into a finite amplitude density hole [65J arising at the scale-size of 
the order of the electron skin depth c/wpe)- 

A powerful circularly polarized electromagnetic (CPEM) plane wave interacting nonlinearly with the EPOs would 
generate an envelope of the CPEM vector potential Aj^ — A±{x+iy) exp{—iLUQt+ikQz), which shall obey the nonlinear 
Schrodinger equation |178j 



2ino 



d_ 
di 



A^ 



7 



where the electron wave function ^ and the scalar potential are governed by, respectively, 

9V Hi d^i^ 



dt 



and 



2 dz^ 



+ (</.- 7 1)V^ = 0, 



- 1 M_L = 0, (130) 



(131) 



IVp-1, (132) 

term. Here fip represents the CPEM wave 
frequency, Vg is the x component of the group velocity of the CPEM wave, iJe is a quantum coupling parameter, and 



where the electron number density is defined as is given by the j'i/'P 

7 = (1 + |Aj^p)^/2 is the relativistic gamma factor due to the electron quiver velocity in the CPEM wave fields. The 
details of normalization of variables is given in Ref. [65j . The nonlinear coupling between intense CPEM w aves and 
EPOs comes about due to the nonlinear current density, which is represented by the term \ip\'^A±/^ in Eq. (130). In 
Eq. (1311, 1 — 7 is the relativistic pondcromotivc potential |178j . which arises due to the cross-coupling between the 
CPEM wave-induced electron quiver velocity and the CPEM wave magnetic field. 

The effect of quantum dispersion on localized electromagnetic pulses can be studied by considering a steady state 
structure moving w ith a constant speed Vg. Inserting the ansatz A±^ — W{£) exp(— ifif), -ip — P{£,) exp{ikx — iuit) and 
(j) = 0(^) into Eqs (fl30| ) -([T32| ), where ^ = z - Vgt, k = Vg/He and uj = /2H^, and where W{i) and P{C) are real, 
one obtains from (ISORlferthe coupled system of equations 



(fW_ 



X- 



pi 



1 W^ = 0, 



(133) 



where 7 = (1 -|- W'^Y^'^, and 



2 de 



7+1)^ = 0, 



(134) 



^P' 



(135) 



with the boundary conditions W" = <i> = and P^ 
frequency shift of the CPEM wave. In the limit Hg 



1 at 1^1 = 00. In Eq. (1331, A = 2r2of^ represents a nonlinear 
0, one has from ( 134 1 (f) — where P ^ 0, and one recovers 

the classical (non-quantum) case of the relativistic solitary waves in a cold plasma |217j . The system of equations 
( 133 1-( 135 1 admits a Hamiltonian 



1 fdW 

2 ['W 



Hi 
2 



dP 



d(j)_ 



1 



(A + l)W'^ + P^ - 7P2 + _ ^ ^ 



(136) 



where the boundary conditions d/d^ — 0, W — (/) — and |P| = 1 at |^| = 00 have been used. 

Numerical solutions of the quasi-stationary system ( 133 1 -(|135 ) are presented Figs. 17 and 18 while time-dependent 



solutions of Eqs. (130l-(132l are shown in Figs. 19 and 1201 Here parameters were used that are representative of 



the next generation laser-based plasma compression (LBPC) schemes [T31 [TH |5T]. The formula |178j eA±/mc'^ — 
6 X 10~^"As"\// will determine the normalized vector potential, provided that the CPEM wavelength As (in microns) 
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and intensity / (in W/cm^) are known. It is expected that in LBPC schemes, the electron number density no may 
reach 10^^ cm~^ and beyond, and the peak values of eA±/m(? may be in the range 1-2 (e.g. for focused EM pulses 
with A, - 0.15 nm and / - 5 X 10^'^ W/cm^). For uOp^ = 1.76 x 10^* s 



one has TiuOpe = 1.76 x 10 erg and 



1 X 10 ^ erg. The electron skin depth Ae ^ 1.7 A. On the other hand, a higher value of 
5.64 X IQ^^ s~^. Thus, our numerical solutions below, based on these two values of 



He = 0.002, since mc^ = i 
Hf, = 0.007 is achieved for = 5.b4 x iU"" s 
He, have focused on scenarios that are relevant for the next generation intense laser-solid density plasma interaction 
experiments [HI]. 

Figures 17 and 18 show numerical solutions of Eqs. (133l-(135) for several values of iJg. The nonlinear boundary 
value problem was solved with the boundary conditions W — (f) — Q and P = 1 at the boundaries at ^ = ±10. One 
can see that the solitary envelope pulse is composed of a single maximum of the localized vector potential W and a 
local depletion of the electron density P^, and a localized positive potential at the center of the solitary pulse. The 
latter has a continuous spectrum in A, where larger values of negative A are associated with larger amplitude solitary 
EM pulses. At the center of the solitary EM pulse, the electron density is partially depleted, as in panels a) of Fig. 
[TT] and for larger amplitudes of the EM waves one has a stronger depletion of the electron density, as shown in panels 
b) and c) of Fig. 17 For cases where the electron density goes to almost zero in the classical case [217], one important 
quantum effect is that the electrons can tunnel into the depleted region. This is seen in Fig. [18) where the electron 
density remains nonzero for the larger value of H^ in panels a) , while the density shrinks to zero for the smaller value 
of He in panel b). 



a) X= 



b) I. 



c) \^-0A 




FIG. 17: The profiles of the CPEM vector potential W (top row), the electron number density (middle row) and the 
scalar potential <1> (bottom row) for A = —0.3 (left column), A = —0.34 (middle column) and A = —0.4 (right column), with 
He = 0.002. After Ref. 65 . 



diffraction effects on the dynamics of localized CPEM wavepackets. 



Figures 19 and 20 show numerical simulation results of Eqs. (130)-(132l, in order to investigate the quantum 

Here the long-wavelength limit o^o ~ 1 and 
Aq — \ and a uniform 



« was considered. In the initial conditions, an EM pump with a constant amplitude A± 
plasma density ip = 1 was used, together with a small amplitude noise (random numbers) of order 10^^ added to 
A± to give a seeding any instability. The numerical results are displayed in Figs. 19 and [20] for He = 0.002 and 
He = 0.007, respectively. In both cases, one can see an initial linear growth phase and a wave collapse at i ~ 70, 
in which almost all the CPEM wave energy is contracted into a few well separated localized CPEM wave pipes. 
These are characterized by a large bell-shaped amplitude of the CPEM wave, an almost complete depletion of the 
electron number density at the center of the CPEM wavepacket, and a large-amplitude positive electrostatic potential. 
Comparing Fig. [l9]with Fig. |20| one can see that there is a more complex dynamics in the interaction between the 
CPEM wavepackets for the larger He = 0.007, shown in Fig. [20] in comparison with He — 0.002, shown in Fig. [19] 
where the wavepackets are almost stationary when they are fully developed. 



IX. SUMMARY AND PERSPECTIVES 



In this paper, we have presented our up-to-date theoretical knowledge of nonlinear physics of non-relativistic 
quantum plasmas. We started with nonlinear quantum models that describe the physics of localized excitations in 
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FIG. 18: The profiles of tlie CPEM vector potential W (top row), the electron number density (middle row) and the scalar 
potential $ (bottom row) for He = 0.007 (left column) and = 0.002 (right column), with A = -0.34. After Ref. gSj. 




20 40 60 



FIG. 19: The dynamics of the CPEM vector potential A± and the electron number density 1-01^ (upper panels) and of the 
electrostatic potential $ (lower panel) for = 0.002. After Ref. |65) . 



different areas of physics. We then moved to discussing the well known electron and ion plasma wave spectra in 
dense quantum plasmas, and presented nonlinear models for treating nonlinear interactions among finite amplitude 
plasma waves at nanoscales. As examples, we demonstrated the existence of localized nonlinear EPOs and localized 
ion waves in dense quantum plasmas. For electrostatic EPOs, the electron dynamics is governed by a pair composed 
of nonlinear Schrodinger and Poisson (NLSP) equations. We stress that the nonlinear Schrodinger equation governing 
the spatio-temporal evolution of the electron wave function in the presence of self-consistent electrostatic potential 
in dense plasmas has been obtained from the quantum electron momentum equation by introducing the eikonal 
representation as a mathematical tool for understanding the complex electron plasma wave interactions at nanoscales. 
Such a description also appears in the context of the nonlinear electron dynamics in thin metal films |112j . Our 
NLSP equations admit a set of conserved quantities, e.g. the total number of electrons, the electron momentum, 
the electron angular momentum, and the electron energy. We have found that the NLSP equations admit quasi- 
stationary, localized structures in the form of one-dimensional quantized dark solitons and two-dimensional quantized 
vortices. These nanostructures are associated with a local depletion of the electron density associated with positive 
electrostatic potential, and are parameterised by the quantum coupling parameter only. In the two-dimensional 
geometry, there exist a class of vortices of different excited states (charge states) associated with a complete depletion 
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20 10 



FIG. 20: The dynamics of the CPEM vector potential A± and the electron number density \tp\^ (upper panels) and the 
electrostatic potential (f> (lower panel) for = 0.007. After Ref. [65] . 



of the electron density and an associated positive potential. Numerical simulations of the time-dependent NLSP 
equations demonstrate the stability of stable dark solitons in one space dimension with an amplitude consistent with 
the one found from the time-independent solutions. In two-space-dimensions, the dark solitons of the first excited 
state were found to be stable and the preferred nonlinear state was in the form of vortex pairs of different polarities. 
One-dimensional dark solitons and singly charge two-dimensional vortices are thus long-lived nonlinear structures in 
quantum plasmas. 

Similar to the Bernstein-Green-Kruskal modes |;164J , we accounted for the trapping of electrons in the electrostatic 
wave potential and studied numerically the deformation of the equilibrium Fermi-Dirac distribution function and the 
subsequent emergence of localized phase space kinetic structures. Such an investigation was carried out by using 
the time-dependent Wigner and Poisson equations. Furthermore, we presented two mechanisms for the generation 
of magnetic fields in quantum plasmas. They are associated with the quantum Weibel instability in the presence of 
equilibrium anisotropic Fermi-Dirac electron distribution function, as well as with the the non-stationary ponderomo- 
tive force of a large amplitude electromagnetic waves in quantum plasmas with streaming electrons. Spontaneously 
generated magnetic fields can affect the linear and nonlinear propagation of both the electrostatic and electromagnetic 
waves in quantum magnetoplasmas [35]. The quantum corrections produce dispersion at short scales for the elec- 
trostatic upper-hybrid, lower- hybrid, and ion-cyclotron waves, while the quantum Bohm force and the electron— 1/2 
spin effects introduce new features to the elliptically polarized extraordinary electromagnetic mode |40j . Furthermore, 
the electron-l|2 spin is responsible for a new perpendicularly propagating (with respect to the magnetic field direc- 
tion) electromagnetic spin wave, which can be excited by intense neutrino bursts in supernovae. The newly derived 
nonlinear EMHD and Hall-MHD equations can be used for investigating the multi-dimensional linear and nonlinear 
electromagnetic waves in quantum plasmas. Finally, we have also presented theoretical and numerical studies of 
stimulated Raman and Brillouin scattering instabilities of a large amplitude electromagnetic wave, and the trapping 
of arbitrary large amplitude circularly polarized EM waves in a fully nonlinear electron density hole in an unmagne- 
tized quantum plasma. It is expected that localized nanostructures can transport electromagnetic wave energy over 
nanoscales in laboratory and astrophysical dense plasmas with degenerate electrons. 

The field of the nonlinear quantum plasma physics is extremely rich and vibrant today, and it holds a great promise 
of providing new practical technologies. For example, the plasma assisted carbon nanostructures and nanomaterials 
are the future of nanotechnologies, so are the new radiation sources in the x-rays and gamma-rays regimes. In such 
circumstances, one should have a complete understanding of the fundamentals of collective nonlinear interactions 
(e.g. intense high-order harmonic generation of ultrashort laser pulses from laser-irradiated dense plasma surfaces) 
in quantum plasmas. Furthermore, in magnetars, and in the next generation intense laser-solid density plasma 
interaction experiments, one would certainly have degenerate positrons, besides degenerate electrons. The physics of 
dense quantum magnetoplasmas with degenerate electron-positron pairs is expected to be quite different than what 
has been described in this paper. The reason is the complex nonlinear dynamics of the electron-positron pairs, which 
would have relativistic velocities in a dense magnetoplasma. Accordingly, we should develop new theories involving 
relativistic kinetic and quantum relativistic magnetohydrodynamic equations that include quantum relativistic effects 
[2181 12191 1220] ■ electromagnetic forces, angular momentum spin and nonlinear effects on equal footings. A detailed 
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analysis of such theories would provide us a guide line for understanding the origin of localized high-energy radiation 
and other complex phenomena (e.g. the formation of structures) from astrophysical settings and future laboratory 
experiments aiming to model astrophysical scenarios. 
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APPENDIX A: DERIVATION OF THE VLASOV EQUATION FROM THE WIGNER EQUATION 

We here show that the Wigner equation converges to the Vlasov equation in the classical limit fi ^ 0. The Wigner 
equation reads 



9/ 
dt 



VV/: 



j_(v- V ) • A 



By the change of variables A = h^/rrie we have d'^A = TC'dPS^/m^ and Eq. ( Al I takes the form 
df 
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Assuming that h/rrie is small and Taylor expanding up to the third order around x, 
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exp[i(v - v') • = iVv' exp[i(v - v') • 
where Vv' = x9/5u^ + yd/dVy + zd/dv'^, we have 
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where the arrows indicate the direction of operation of the nabla operators. 

Integration by parts in v' space, and using that / — > for |v'| oo, now yields 



df 



vV/ = -- 



d-^^d^v' exp[i(v - v') • $,] { 0(x, t) 



dt ' (27r)3m 

Now the integration in ^ space can formally be performed, with the result 
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was used and where S is Dirac's delta function. Finally, integration over v' space yields 
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In the limit ?i — s- 0, we recover the Vlasov equation 
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APPENDIX B: DERIVATION OF THE DISPERSION RELATION FOR THE WIGNER-POISSON 

SYSTEM 

We here present a derivation of the dispersion relation for electrostatic waves in a degenerate quantum plasma, 
governed by the Wigner-Poisson system of equations. The linearized Wigner-Poisson system of equations reads 



dt 



V V/i = 



(27r)3?i4 



and 
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where /o denotes the background distribution function, and /i and (pi denote the perturbed distribution function and 
electrostatic potential, respectively. Fourier transforming (B1|-(B2) in space and Laplace transforming them in time, 
we obtain 
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Rewriting ( B3 1 as 
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and performing the integration over A space, we have 
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where 5 is the Dirac delta function. Now, the integration can be performed over v' space, obtaining the result 
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Solving for /i in ( B7 1 and inserting the result in (|B4|) , we obtain the dispersion relation 
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Suitable changes of variables in the two terms now gives 
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which can be rewritten as 



Aire'' 



/o(u) 



(^-k.u)2-|^ 



d^u = 0. 



(B6) 



(B7) 



(B8) 



(B9) 



(BIO) 



The dispersion relation (BIO) was also derived by Bohm and Pines [3J using a series of canonical transformations of 

n (57) 
^hen, ^ 

47re^ f /o(u) 



the Hamiltonian of the system; see for example equation (57) in their paper. We now choose a coordinate system 
such that the x axis is aligned with the wave vector k. Then, (BIO) takes the form 
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(Bll) 
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We next consider a dense plasma with degenerate electrons in the zero temperature limit. Then, the background 
distribution function takes the simple form 



1 0, elsewhere, 



(B12) 



where Vp^ — (2£i?e/me)^/^ is the speed of an electron on the Fermi surface, and Spe = (37r^no)^/'^?i^/(2TOe) is the 
Fermi energy. The integration in (Bill can be performed over velocity space perpendicular to u^, using cylindrical 
coordinate in Uy and u^, obtaining the result 
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Equation (B13l can be written as 
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First, in the limit fik/m^ 0, we have from (B15l 
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where we have assu med that lo is real and io/k > Vpe- 

Expanding (B16l for small wavenumbers up to terms containing fc^, we have 
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Second, for non-zero hk /m^, we have from (B15) the dispersion relation 
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(B18) 



Expanding (B18l for small wavenumbers up to terms containing fc**, we obtain 



(B19) 



where a = {48/175)mlVk/n^ojle ~ 2mO{nQaiy/^ and = r? jmef?- « 53 x 10"^° cm is the Bohr radius. For a 
typical metal such as gold, which has a free electron number density of no = 5.9 x 10^^ cm^'^, we would have a k, 0.4. 
For the free electron density in semiconductors, which is many orders of magnitude less than in metals, a is much 
smaller and can safely be dropped compared to unity. 
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